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Abstract 

An objective of the High Performance Computing 
and Communication Program at the NASA Langley 
Research Center is to demonstrate multidisciplinary 
shape and sizing optimization of a complete aerospace 
vehicle configuration by using high-fidelity, finite- 
element structural analysis and computational fluid 
dynamics aerodynamic analysis. In a previous study, a 
multi-disciplinary analysis system for a high-speed civil 
transport was formulated to integrate a set of existing 
discipline analysis codes, some of them 
computationally intensive. This paper is an extension 
of the previous study, in which the sensitivity analysis 
for the coupled aerodynamic and structural analysis 
problem is formulated and implemented. Uncoupled 
stress sensitivities computed with a constant load vector 
in a commercial finite element analysis code are 
compared to coupled aeroelastic sensitivities computed 
by finite differences. The computational expense of 
these sensitivity calculation methods is discussed. 

Nomenclature 

a c Aerodynamic cruise load vector 

a M Aerodynamic maneuver (converged) load vector 

f A Augmented structural load vector 

f c Structural cruise load vector 

f M Structural maneuver (converged) load vector 

s 0 Unloaded shape vector 

s c Cruise shape vector 

s M Maneuver (converged) shape vector 

<t a Stress vector (due to augmented loads) 

u 0 Displacement (unloaded to cruise shape) vector 

u A Augmented displacement vector 

u c Cruise displacement vector 

ii M Maneuver (converged) displacement vector 

v Independent design variables 


w c Cruise weight 
w M Maneuver weight 
K c Stiffness matrix (cruise shape) 

K 0 Stiffness matrix (unloaded shape) 

LF Inertial load factor (g’s) 

Introduction 

One of the objectives of the High Performance 
Computing and Communication Program (HPCCP) at 
NASA Langley Research Center (LaRC) has been to 
promote the use of advanced computing techniques to 
rapidly solve the problem of multidisciplinary 
optimization of aerospace vehicles. In 1992, the 
HPCCP Computational Aerosciences (CAS) team at 
LaRC began a multidisciplinary analysis and 
optimization software project. Initially, the focus of the 
CAS project was on the software integration system 
used to integrate fast analysis on a simplified design 
application. The sample application for this project was 
a High Speed Civil Transport (HSCT). Over the years, 
progressively more complex engineering analyses have 
been incorporated. In 1997, the sample application 
shifted to a more realistic model and higher fidelity 
analyses and is referred to as the HSCT4.0 application. 
The analysis formulation and preliminary results from 
the HSCT4.0 application were presented in Refs. 1 and 
2, respectively. However, sensitivity analysis was not 
performed in Refs. 1 and 2. 

As discussed in Refs. 1 and 2, analytical 
sensitivities can be obtained from most of the analysis 
codes by using automatic differentiation tools. 
However, Ref. 1 discussed the one major stumbling 
block in formulating the sensitivity analysis — obtaining 
sensitivity derivatives of the converged aeroelastic 
loads and using those derivatives in a commercial finite 
element code. The aerodynamic pressures and 
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structural deformations in the wing shape are mutually 
dependent. Because of this coupling between 
aerodynamic and structural responses, the sensitivity of 
the aeroelastic loads is difficult to compute. 
Researchers use one of two methods to compute these 
aeroelastic sensitivity derivatives. For the first method, 
researchers assume that the aeroelastic loads are 
constant and then use exiting methods to obtain 
sensitivity derivatives. In the second method, 
researchers assume the loads vary and use finite 
differences to obtain sensitivity derivatives. There has 
been some research 2 3 * 5 using the global sensitivity 
equations (GSE 6 ) to account for the coupling of 
aerodynamics and structures. In Refs. 4 and 5, the GSE 
were applied to simplified HSCT models to compute 
derivatives of aerodynamic coefficients. In this paper, 
the effects of the coupled aeroelastic loads on the 
computation of stress sensitivities will be examined. 

In Refs. 1 and 2, the HSCT4.0 analysis consisted of 
an integrated set of discipline codes and interface 
codes. These codes were implemented in a CORBA- 
Java computing environment known as CJOPT 7 . The 
present research does not use the CJOPT system; 
instead, the discipline codes and interface codes are 
implemented in the Phoenix Integration, Inc/s 
ModelCenter® 8 . Implementation of the codes in 
ModelCenter® is not discussed in this paper. The work 
presented in the current paper is referred to as the 
HSCT4.1 application. 

In the HSCT4.1 application, the coupling of the 
aeroelastic loads in the stress derivatives is important 
because both shape and structural design variables are 
used in the sensitivity analysis. The present paper will 
quantify the effect of aeroelastic coupling on stresses. 
Coupled stress sensitivities obtained by finite difference 
techniques (with loads varying) are compared with un- 
coupled stress sensitivities obtained from a commercial 
finite element code (with loads held constant). 

First, an overview of the model and analysis will 
be presented. Next, the coupled and uncoupled stress 
sensitivity analysis formulation will be presented. 
Finally, sample results for both the coupled and 
uncoupled stress sensitivities are compared for both 
shape and structural design variables. 

Overview 

HSCT 4.1 Model 

The HSCT4.1 application uses the same linear 
aerodynamic and finite element models as the HSCT4.0 
application described in Refs. 1 and 2. The linear 
aerodynamics grid and the finite element model used 
for HSCT4.1 are shown in Figs. 1 and 2, respectively. 
An aerodynamic surface grid of approximately 1100 
grid points is used in the linear aerodynamics code 
(USSAERO) 9 . The aerodynamic model does not 


include engines. The GENESIS 010 finite element 
analysis code (a product of VR&D, Inc.) uses the 
40,000 degree- of- freedom (DOF) FEM for 

displacement and stress response calculations. In this 
FEM, the engines are modeled as masses on beam 
elements. Seven laterally symmetric load conditions 
are used — one cruise load condition and six maneuver 
conditions (three at +2.5g and three at -i.Og). Note that 
the taxi condition used in Refs. 1 and 2 is not used in 
the present work. 



Fig, 1. HSCT4,1 linear aerodynamic model. 



Fig. 2. HSCT4.1 Finite element model. 


Optimization Problem Description 

The objective function of the HSCT4.1 
optimization problem is to minimize the aircraft gross 
takeoff weight (GTOW) subject to stress constraints. 
The HSCT4.1 application has 271 design variables for 
optimization — 244 structural thickness variables and 27 
shape variables. To limit the number of independent 
structural design variables, the optimization model is 
divided into 61 design variable zones. Each zone 
consists of several finite elements. Thirty-nine zones 
are located on the fuselage, and twenty-two zones are 
located on the wing (half are on the upper surface and 
the other half are on the lower surface). Within each 
zone, four structural design variables are used. These 
structural design variables consist of three ply-thickness 
variables (a 0° fiber variable, a 45° fiber variable, and a 
90° fiber variable with a dependent -45° fiber variable 
set equal to the 45° variable) and a core thickness 
variable. The 0° and 90° ply orientations at various 
locations on the model and the composite laminate 
stacking sequences are shown in Fig. 3. 

The 27 shape design variables are divided into two 
sets. The first set contains 9 planform design variables- 
the root chord, the outboard break chord, the tip chord, 
the distance from the semispan to the outboard break, 
the leading edge sweep of the inboard wing panel, the 
leading edge sweep of the outboard wing panel, the 
total projected area of the wing, the fuselage nose 
length, and the fuselage tail length. The second set 
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consists of three sets of six variables that control the 
camber, thickness, and shear of the wing airfoil. 




Fig. 3. Ply orientations and composite laminate 
stacking sequence. 


HSCT4. 1 Analysis Process 
As discussed in Refs. 1 and 2, the HSCT4.0 
analysis is formulated as a sequence of processes in the 
data flow diagrams, as shown in Fig. 4. The HSCT4.I 
analysis is focused on the difficult aerodynamic- 
structures analysis-coupling problem and uses the 
analysis processes shown by shaded circles in Fig. 4. 
By convention, this paper uses italics for process names 
from Fig. 4. 



Fig. 4. Analysis Process. 


The HSCT 4.1 Analysis process is started from the 
top of the data flow (Fig. 4) where the design variable 
values are prescribed. First, the Geometry process is 
used to derive updated geometric grids and FEM 
section properties from the design variables. Next, the 
Weights process uses the derived FEM grid and section 
properties to calculate detailed weights and the center 
of gravity (c.g.) locations for specified flight conditions. 
Theoretical FEM weights are computed for each node 
from the FEM data. A reference as-built weight 
increment is added to the theoretical FEM weight at 
each node to produce the as-built weight. 

Next the Rigid Trim process is executed for the 
cruise condition to determine the configuration angle of 
attack and the tail deflection angle that combine to yield 
a lift equal to the weight, with no net pitching moment. 
The linear aerodynamics code is used to compute the 
trimmed aerodynamic pressures on the cruise shape 
aerodynamic grid. These pressures are then transferred 
from the aerodynamic grid points to the FEM nodal 
forces in the z-direction. Inertial forces (nodal weights 
times g-factor) are added to the aerodynamic forces to 
create the structural load vector. This structural load 
vector is used by the Displacements process to compute 
the unloaded shape of the aircraft (see Appendix A for 
further details). These cruise displacements are saved 
as a reference set for the Loads Convergence process. 

In the Loads Convergence process, aeroelastic trim 
calculations are performed for the six noncruise load 
conditions, to produce the aeroelastically converged 
loads on the aircraft (see Appendix B for further 
details). Forces representing cabin pressure are added 
to these converged loads and are multiplied by a factor 
of safety (1.5). These augmented loads are used in the 
Stress & Buckling process to compute stress and buck- 
ling constraints for all elements contained in the 61 
design zones on the fuselage and wing. Buckling 
constraints are not considered in this paper. 

HSCT4.1 Sensitivity Analysis Formulation 

In this section, sensitivities (first order derivatives 
with respect to the design variables) are formulated for 
the objective function (GTOW) and the stress con- 
straints. First, the weight sensitivities are formulated 
analytically. Then the stress sensitivities for both the 
coupled and uncoupled methods are formulated. 

Weight Sensitivity Formulation 

In the HSCT4.1 system, structural weight sensi- 
tivities are computed analytically from the finite 
element dimensions and material properties using chain 
rule differentiation. Non-structural weights (fuel, 
payload, and non-modeled structures) are assumed to be 
constant (Le. non-structural sensitivities are zero). 

Only shell elements are sized by the structural 
design variables. Because the weight of each of these 
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elements is a linear function of the structural design 
variables, the weight sensitivity for those elements with 
respect to a structural design variable is simply the 
element area multiplied by the number of plies and the 
density of that ply. 

Element weight sensitivities with respect to the 
shape design variables are computed by chain rule dif- 
ferentiation of the equations defining the dimensions of 
the element (length for beams, area for shells, and vol- 
ume for solid elements). The Geometry process is used 
to compute the derivatives of the node locations with 
respect to the shape design variables. Element weights 
and weight sensitivities are extrapolated to the finite 
element nodes by the weight code. The element weight 
sensitivities are totaled to obtain the GTOW sensitivity. 


Stress Sensitivity Formulation 

Structural stresses, C A , are functions of the 
displacements. Therefore, the stresses and stress 
sensitivity derivatives (da A ldv) are defined as follows: 

°A = o A (u A ) ( 1 ) 


d<y A 

dv 


dga 

du A dv 


( 2 ) 


The relationship between the stresses and the displace- 
ments are defined in GENESIS®; therefore, it is only 
necessary to study the displacement derivatives 
(dujdv). Linear static structural displacements at 
converged maneuver conditions ( u A ) are computed 
within GENESIS® from the stiffness matrix (K c ) and 
the augmented loads (f A ) using Equations 3 and 4. 

= fA (3) 

= [KcV/a W 

Derivatives of u A with respect to the design variable 
vector (v) are obtained by differentiating Equation 3: 


K/ U * 


dv 


dKc 

dv 


dv 


(5) 


Rearranging Equation 5 gives: 


du^ 

dv 


= [K C V 


dh 

dv 


[K c Y dKc 


dv 


(6) 


In Equation 6, the displacement derivative is composed 
of two terms. One term requires the force derivative, 
which is derived from the coupled aerodynamic and 
structural analyses, and is referred to as the coupled 
sensitivity term in Equation 7. 


d“ A 

dv 


JCoupied 



(7) 


The other term requires the stiffness matrix derivative, 
which is determined from structural terms only, and is 
referred to as the uncoupled sensitivity term, Equation 
8 . 


^ = - [K C V^U A ( 8 ) 

dv „ , , dv 

J Uncoupled 

The relationship between structural design vari- 
ables and the stiffness matrix can be set up as direct 
input to GENESIS®. The Geometry > process was used 
to generate relationships between the nodes of the stiff- 
ness matrix and the shape design variables as another 
set of GENESIS® input. Therefore, the uncoupled 
sensitivity term (Equation 8) can be computed directly 
using the sensitivity analysis capability in GENESIS^. 

The coupled sensitivity term (Equation 7) is similar 
to the displacement calculation equation (Equation 4). 
Therefore, GENESIS® can be used to compute the 
coupled sensitivity term in the displacement sensitivity 
equation from a linear static finite element displace- 
ment analysis using the load sensitivities as a nodal 
force vector. One complication of this method is that 
the sensitivities with respect to each design variable 
must be defined as a separate load case (i.e. one load 
case is required per design variable). Another compli- 
cation is the fact that the augmented load sensitivity 
(df A fdv) must be computed outside of GENESIS®. 

The augmented loads (f A ) are computed from the 
converged maneuver loads (f M ), the factor of safety 
(1.5) and the cabin pressures (constant); therefore the 
derivatives of the augmented loads are: 

4f A _ 1 r M (Q\ 

dv dv 

In the Loads Convergence process, the maneuver 
loads (f M ) and maneuver displacements (u M ) are 
mutually dependent. The nature of this dependency is 
discussed in Appendix B, and is represented in 
Equations 10 and II: 

= «m(/ M ’ V ) ( 10 ) 

f M = ( U ) 

Because of the coupling of f M and u M , computation 
of df M /dv is complicated and beyond the scope of this 
paper. The total stress sensitivity derivative (da A ldv) is 
obtained by substituting du A ldv from Equation 6 into 
Equation 2. For the uncoupled stress sensitivity 
derivative, du A ldv from Equation 8 are used in Equation 
2 instead. In the next section, results are presented for 
analytically computed uncoupled stress sensitivities. 

Coupling effects are shown by parametric stress 

analyses in the next section. 

HSCT4.1 Weight and Stress Results 

Analysis results were presented in Ref. 2. To 
validate the HSCT4.1 Analysis process, loads 

convergence displacements and stress results from Ref. 
2 for the baseline configuration were compared with 
HSCT4.1 results. The results from the two analyses 
agreed with each other. In Ref. 2, ten iterations were 
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used in the Loads Convergence process. Therefore, all 
results in this paper used ten iterations in the Loads 
Convergence process. Both parametrically- and 
analytically-computed GTOW and stress results are 
presented in this section. 

Parametric studies were performed by perturbing 
one design variable at a time from its baseline value and 
then running the complete HSCT 4.1 Analysis . In the 
parametric studies, design variables are perturbed by a 
maximum of 1% from their baseline values. 

Structural responses (r, referred to as analytical 
responses) are computed from Equation 12. Sensitivity 
analyses are used to compute first order derivatives of 
the response functions ( drldv ) at the baseline design 

point {y Baseline)* 

K v ) = ~ ( V “ V Baseline ) + r Baseline 

dv 

Due to the large number of design variables 
available for study, results for a small subset of the 
design variables are presented in this paper; four of the 
twenty-seven shape design variables (two chord 
lengths, one span length, and the wing planform area) 
are selected for study. Further, only two of the 244 
structural design variables (ply thicknesses in the 
design variable zone at the wing break where the largest 
stresses occurred) are selected for study. 

Weight Variations with Design Variables 

Weight results are presented as the GTOW 
normalized by the baseline GTOW. GTOW is plotted 
as a function of the six selected design variables in 
Figs. 5 to 10. The analytic GTOW curves are computed 
using Equation 12 (linear approximations). Note that 
the weight sensitivity calculations do not involve 
aeroelastic coupling. The analytic curves agreed with 
the parametric curves. The GTOW axis is different in 
Figs. 5 to 8 (variations in shape variables) due to the 
large difference in the GTOW variations for each 
variable. As described in the Weight Sensitivity section 
and shown in Figs. 9 and 10 (variations in structural 
variables), weight is a linear function of the each of the 
structural design variables. 



Fig. 5. GTOW variation with design variable 3. 



Fig. 6. GTOW variation with design variable 4. 



Fig. 7. GTOW variation with design variable 5. 
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Fig. 8. GTOW variation with design variable 7. 
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Fig. 9. GTOW variation with design variable 184. 



Design Variable 185 (Normalized) 

Fig. 10. GTOW variation with design variable 185. 

Stress Variations with Design Variables 

Stress results are presented as stress failure index 
(SFI) values (Ref. 2). In the Analysis process, stresses 
are computed for eight plies of 2260 sized elements for 
six load conditions. Due to the large number of results 
(108,480 responses), only subsets of these results are 
presented in this paper. Stress results are shown only 
for one of the three +2.5g maneuver load conditions 
and for three elements (see Fig. 11) with large stress 
responses. 



Fig. 11. Elements selected for stress responses. 


In Figs. 12 to 29, parametrically and analytically 
computed SFI values are plotted as functions of the 
design variables. All analytical SFI curves are 
determined using Equation 12 where the uncoupled 
stress sensitivity (obtained from GENESIS 0 using 
Equations 2 and 8) is used for the response sensitivity 
term ( drfdv ). In the following plots this stress 

sensitivity is represented as the slope of the analytic 
stress curve. 

SFI responses for element 5140 (on the fuselage) 
are plotted in Figs. 12 to 17. The ranges of the SFI axes 
in Figs. 12 to 15 are identical for comparison purposes. 
The parametric SFI response curves in Figs. 12 to 15 
are shallow curves (not quite linear). It is noted that the 
slope of the parametric curves at the baseline design 
point (normalized design variable value of 1.0) are 
different from the slope of the analytic (uncoupled) 
curve. The difference in the slopes varies from one 
variable to another, but in general the difference is 
significant. The ranges of the SFI axes in Figs. 16 and 
17 are magnified from Figs. 12 to 15 because the 
sensitivity of the SFI in the fuselage is much smaller 
with respect to changes in the upper wing surface (zone 
40) ply thicknesses. Because the sensitivity with 
respect to variables 184 and 185 is very small, the SFI 
response is more erratic. Some of this noise is the 
result of round-off error due to passing results from one 
analysis code to another using low-precision text 
output. This noise could also be due to the nonlinearity 
of the SFI response. 

Note that the differences in the slopes of the 
parametric and analytic (uncoupled) SFI response 
curves are larger in Figs. 12 to 15 than in Figs. 16 and 
17. Thus, the coupling effect (due to aeroelastic forces) 
seems to be larger for the shape design variables than 
for the structural design variables. The aeroelastic 
forces are a function of the weight and the HSCT shape. 
Shape variables affect both the shape and the weight, 
while the structural design variables only affect the 
weight. Therefore, the shape variables have a larger 
coupling effect for the SFI responses than the structural 
variables; a conclusion, which is supported by Figs. 12 
to 17. 
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Fig. 12. SFI variation with design variable 3 for 
element 5140. 
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Fig. 15. SFI variation with design variable 7 for 
element 5140. 
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Fig. 13. SFI variation with design variable 4 for 
element 5140. 
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Design Variable 184 (Normalized) 

Fig. 16. SFI variation with design variable 184 for 
element 5140. 



Fig. 14. SFI variation with design variable 5 for 
element 5140. 



0.990 0.995 1.000 1.005 1.010 

Design Variable 185 (Normalized) 

Fig. 17. SFI variation with design variable 185 for 
element 5140. 


SFI responses for element 25035 (on the upper 
wing surface, see Fig. 11) are plotted in Figs. 18 to 23. 
The ranges of the SFI axes in Figs. 18 to 21 are 
identical for comparison purposes. The parametric SFI 
responses in Figs. 18 to 21 are nearly linear functions 
within the selected design variable range. Again, the 
slopes of the parametric SFI curves are different from 
the slopes of the analytic (uncoupled) SFI responses. 
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The ranges of the SFI axes in Figs. 22 and 23 are 
different than Figs. 18 to 21 because the sensitivities are 
smaller. Unlike the SFI sensitivities in Figs. 16 and 17, 
the sensitivities in Figs. 22 and 23 are larger than the 
numerical noise in the SFI variation. The larger 
sensitivities in Figs. 22 and 23 are because the thickness 
of element 25035 is sized by structural design variables 
184 and 185. As in the earlier figures, the effect of the 
coupling term in the SFI sensitivities is larger for the 
shape variables than for the structural variables. 



Fig. 18. SFI variation with design variable 3 for 
element 25035. 



Fig. 19. SFI variation with design variable 4 for 
element 25035. 



Fig. 20. SFI variation with design variable 5 for 
element 25035. 



Fig. 21. SFI variation with design variable 7 for 
element 25035. 



Fig. 22. SFI variation with design variable 184 for 
element 25035. 



Design Variable 185 (Normalized) 

Fig. 23. SFI variation with design variable 185 for 
element 25035. 

SFI responses for element 25082 (on the lower 
wing surface, see Fig. 1 1) are plotted in Figs. 24 to 29. 
The ranges of the SFI axes in Figs. 24 to 27 are 
identical for comparison purposes. The parametric SFI 
responses in Figs. 24 to 27 are nearly linear functions 
within the selected design variable range. Again, the 
slopes of the parametric SFI curves are slightly 
different from the slopes of the analytic (uncoupled) 
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SFI responses. The ranges of the SFI axes in Figs. 28 
and 29 are magnified from Figs. 24 to 27 because the 
sensitivity of the SFI in the lower wing surface is much 
smaller with respect to changes in the upper wing 
surface (zone 40) ply thicknesses. As before, because 
the sensitivity with respect to variables 184 and 185 are 
very small, the effects of noise in the SFI response is 
more pronounced. As in the earlier figures, the effect 
of the coupling term (Equation 7) in the SFI 
sensitivities is larger for the shape variables than for the 
structural variables. 



Fig. 24. SFI variation with design variable 3 for 
element 25082. 



Fig. 25. SFI variation with design variable 4 for 
element 25082. 



Fig. 26. SFI variation with design variable 5 for 
element 25082. 



Fig. 27. SFI variation with design variable 7 for 
element 25082. 



Fig. 28. SFI variation with design variable 184 for 
element 25082. 
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Fig. 29. SFI variation with design variable 185 for 
element 25082. 


terms in the stress sensitivities and in parallel execution 
oftheHSCT4.1 system. 
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Table 1. CPU Time Required for Analyses 



CPU Execution Time (hrs.) j 

Process 

Analysis 

Finite 

Difference 

Uncoupled 

Stress 



Sensitivity 

Sensitivity 

Aero. 

6.10 

1659.20 

6.10 

Disp. 

0.37 

99.73 

0.37 

Stress 

0.03 

9.07 

0.03 

d Stress/ ov 

0.00 

0.00 

0.50 

Total 

6.50 

1768.00 

7.00 


Concluding Remarks 
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ModelCenter® framework is presented. Formulations 
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Appendix A - Unloaded Shape Computation 
Both the geometry of the aerodynamic model 
used in the Rigid Trim process and the FEM used in the 
Displacements process are based on the 1-g cruise 
shape ( Sc, see Fig. A-l), which is described by the 
shape design variables. 


t t t 



t t t 


Fig. A-l. Cruise Shape. 


In order to obtain accurate stress responses from a 
structural analysis, the displacements must be computed 
relative to the unloaded shape (s 0 , see Fig. A-2). The 
unloaded shape is the aircraft shape that deforms to the 
cruise shape when cruise loads (f c ) are applied. 


*0 



Fig. A-2. Unloaded Shape. 


The cruise loads are computed from Equation A.l: 

fc = a c - W c (A.l) 

where a c represents the aerodynamic cruise loads, and 
w c represents the inertial loads. The inertial loads 
(nodal cruise weight) are subtracted from a c because 
the weight acts in the negative z direction. 

The shape change (displacements) from the 
unloaded shape to the cruise shape is given by the 
following equation: 

J c -5 0 = «0 = KVfc ( A - 2 ) 

Unfortunately, the stiffness matrix of the 
unloaded shape (K 0 ) is unknown. Assuming that the 
FEM is geometrically linear, the differences between 
the stiffness matrices of the unloaded shape (K 0 ) and 
the cruise shape ( K c ) are negligible. Using this 
assumption, the cruise loads (fc) can be applied to the 
cruise shape stiffness matrix ( K c ) to produce a set of 
cruise displacements (m c ). 

u c - [K c \ f c = f c = w 0 (A. 3) 

Thus, the unloaded shape is given by the cruise shape 
( s c ) minus the cruise displacement vector ( u c )• 

Jq — Sq u.£ — s £ Uq (A. 4) 


Appendix B - Loads Convergence 
In this appendix, load conditions 2 through 7 are 
referred to as maneuver conditions, and the 
corresponding loads are converted to forces acting at 
FEM nodes. At maneuver conditions, the aircraft has 
shape s M and is subjected to maneuver loads f M 
(illustrated in Fig. B-l). 



The maneuver loads are computed from Equation B.l : 

/*=«!#- LF W M ( B1 ) 

where a M represents the aerodynamic maneuver loads; 
LF represents the maneuver load factor, and w M 
represents the nodal weight at maneuver conditions 
(GTOW). The Weights process is used to compute the 
nodal GTOW. The inertial loads (LF multiplied by 
nodal GTOW) are subtracted from a M because the 
weight acts along the negative z direction. 

The maneuver shape (s^) and aerodynamic 
maneuver loads (ajf) are unknowns that are determined 
using the Loads Convergence process. The maneuver 
shape is the net aircraft shape used for the aerodynamic 
analysis and consistent with the deformed structural 
shape produced by the aerodynamic loads. 

The maneuver displacements are produced by 
deforming the unloaded shape (sq) to the maneuver 
shape Gm) as illustrated in Fig. B-2. 



The maneuver displacements (um) are computed using 
equation B.2: 

u m = s m~~ s o ~ M 

The unloaded stiffness matrix (K 0 ) is unknown, but by 
applying the geometric linearity assumption from 
Appendix A, the cruise stiffness matrix (K c ) can be 
used in place of K 0 to calculate the maneuver 
displacements (u M ). 

u m ~ s m ~ s o = I^cT ft* 

In the aerodynamics analysis, the maneuver 
shape ( s M ) must be used to compute the aerodynamic 
maneuver forces (a M ). The maneuver shape is obtained 
by perturbing the existing aerodynamic model, which is 
based on the cruise shape (%)- The displacements 
(referred to as "delta displacements," 5 u M ) used to 
perturb the aerodynamic model are computed from 
equation A.4 (in Appendix A) and equation B.3. 

8 U M ~ s m ” s c ~ ( 5 m ~ s o )” i 5 c ~ 4 ^ 

= «M' W C 
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